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•^ ■ Wc consider a Kcllcr-Scgcl model coupled to the incompressible Navier-Stokes equations 

qq , in spatial dimensions two and three. Wc establish the local existence of regular solutions 

£\i ■ and present some blow-up criteria for both cases that equations of oxygen concentration is 

of parabolic or hyperbolic type. We also prove global existence and decay estimate in time 

under the some smallness conditions of initial data. 
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1 Introduction 



In this paper, we consider mathematical models describing the dynamics of oxygen, swimming 
bacteria, and viscous incompressible fluids in R , with d = 2, 3. Bacteria or microorganisms 
often live in fluid, in which the biology of chemotaxis is intimately related to the surrounding 
physics. Such a model was proposed by Tuval et al.|24] to describe the dynamics of swimming 
bacteria, Bacillus subtilis. We consider the following equations in [23] and set Qt = (0, T] xR 
with d = 2, 3: 



' d t n + u ■ Vn — An = —V • (x(c)n'Vc), 

dtc + u-Vc- fiAc = -k(c)n, in (x, t) G M. d x (0, T], (1.1) 

, d t u + u • Vu — Au + Vp = —nV(f>, V • u = 

X' 

where c(t, x) : Q T ->• R+, n(t, x) : Q T -)• R+, u(t, x) : Q T -> R d and p(t,x) : Q T -> R denote 
the oxygen concentration, cell concentration, fluid velocity, and scalar pressure, respectively. 
The nonnegative function k(c) denotes the oxygen consumption rate, and the nonnegative 
function x{ c ) denotes chemotactic sensitivity. Initial data are given by (no(x),co(x),uo(x)). 
We study both cases that either u = 1 or u = in the equation of oxygen and, for convenience, 
the case fi = 1 of (jl.ip is called parabolic Keller-Segel-Navier-Stokes equations (abbreviated 
to P-KSNS) and the case u = is referred as partially parabolic-hyperbolic Keller-Segel- 
Navier-Stokes equations (abbreviated to PH-KSNS). We will refer the system to the Keller- 
Segel-Stokes equations (abbreviated to PH-KSS) if the convection term u ■ Vu is absent in 

(1-1)3. 

To describe the fluid motions, we use Boussinesq approximation to denote the effect due 
to heavy bacteria. The time-independent function <j) = (f>(x) denotes the potential function 
produced by different physical mechanisms, e.g., the gravitational force or centrifugal force. 
Thus, <j)(x) = axd is one example of gravity force, and <j>(x) = <fi(\x\) — > as |x| — > oo is an 
example of centrifugal force. 



The classical model to describe the motion of cells was suggested by PatlakpjJ] and Keller- 
Segel[13j Q3]. It consists of a system of the dynamics of cell density n = n(t,x) and the 
concentration of chemical attractant substance c = c(t, x) and is given as 

{rit = An — V • (riV^c), 
(1-2) 
act = Ac — tc + n, 

where x is the sensitivity and r~2 represents the activation length. The system (|1.2j) has been 
extensively studied by many authors and we will not try to give list of results here (see e.g. 
PH E21 S3 UM ES] and references therein). 

Our main objective of this paper is to present blow-up criteria of (jl.ip in two or three 
dimensions, unless solutions exist globally in time (see Theorem Q] and Theorem [3] below), and 
to establish global existence of regular solutions and their decay properties, when certain norm 
of initial data is sufficiently small (see Theorem [2] and Theorem U] below) . 

We mention previously known results related to ours. In |16j local existence of solutions 
was shown in three dimensional bounded domains and [8] proved the global-in-time existence 
of the smooth solutions when initial data are close to constant states in M 3 and x(')>^(") 
satisfy certain conditions. More precisely, [8] showed that if initial data ||(no — n^, Co, -ucOH/p 
is sufficiently small, then there exists a unique global solution, provided that 

*'(•)> 0, k'(.)>0, (Ml) <0. (1.3) 

In the absence of the fluid in (II. ip . i.e., u = 0, [22] showed that there exists a unique, global 
and bounded solution if x is sufficiently small, dependent upon ||co||£oo( R d). 

For two dimensional case, in [15] . Liu and Lorz showed the global existence of a weak 
solution in IR 2 under the following conditions on x(') an d k(-): 

X'0)>0, (x(0*0))' > 0, (^|) <0. (1.4) 



In two dimensions, Winkler [26] proved the global existence of regular solutions without 
smallness assumptions on initial data for bounded domains with boundary conditions d u n = 
d v c = u = under the following sign conditions on x{') an d k(-): 

(f)'>o. w . w .)y>«, ($)"<o. (i.5) 

In [T] the authors of the paper established global existence of smooth solutions in M. 2 with no 
smallness of the initial data and a certain conditions, motivated by experimental results in [2] 
and [M], on x(-) and k(-) (compare to (|1.5p ). that is, 

x(c), k(c), x'( c )i k'(c) > 0, and sup|x(c) — fJ>k(c)\ < e for some fj, > 0. (1-6) 

Construction of weak solutions in K 3 was also discussed in pQ with replacement of \x( c ) ~ ^k(c)\ - 
in (|1.6|) . We refer to [U] , [23] and [3] and references therein for the nonlinear diffusion models 
of a porous medium type An m , instead of An. 

As mentioned earlier, our main motivation is to study existence of regular solutions of (II. ip 
when certain norm of initial data is small. To be more precise, we show that in case \i = 1, if 



||co||x,oo is small, then solutions become regular in R , d = 2,3 and satisfy a certain degree of 
decay in time (when d = 3, Stokes system is under our consideration for fluid equations). On 
the other hand, in case \x = 0, we establish local solutions in time, and then global solutions 
with time decay if ||no|| d is small. We first consider the case (P-KSF) in (jl.ip . Local-in- 
time existence of classical solutions for (jl.ip was established in Theorem 1 in [T] such that 
(n,c,u) E L oo (0,T;H m ~ 1 (R d ) x H m (R d ) x H m (R d )) for some T > and m > 3. We present 
some blow-up criteria of local classical solutions, unless the maximal existence time is infinite. 

Theorem 1 Let the initial data (n ,c ,ito) be given in H m ~ 1 (R d ) x H m (R d ) x H m (R d ) for 
m > 3 and d = 2,3. Assume that X)k,x'>k' are all non-negative and x, k E C m (R + ) and 
k(0) = ; ||V 0||i°° < oo for 1 < |/| < m. If the maximal time of existence, T* , in Theorem 1 
in fry, is finite, then one of the following is true in each case ofR 2 or R 3 , respectively: 

2 2 
(2D) IMIl 9 (0,T*;Lp(R2)) = °°. - + - = 2, 1 < / < oo. (1.7) 



where 



(3-D) II u IIlt(0,T*;L^(K 3 )) + \\ n \\L9(0,T*;LP(U 3 )) ~ °°' 0-°) 

3 2 3 2 3 

-5 + -<l, 3</3<oo, - + - = 2, -<p<oo. 
P 7 P 9 2 

If fluid equation is the Stokes system for (3D), \\u\\ L gr 0T „. LP , R:i ^ in (jl.8p is dropped. 

The proof of Theorem [1] will be given in section 2. 

Remark 1 VKe remind the following scaling invariance of (11. ip .- 

n R (£, x) := R 2 n(R 2 t, Rx), c R (t, x) = c(R 2 t, Rx), u R (t, x) = Ru(R 2 t, Rx) (1.9) 

and observe that (jl.7p and (|1.8j) are invariant functionals under the scaling (jl.9j) . For t/ie limit- 
ing case (l,m) = (1, oo) in (|1.7|) . due to conservation of total mass, we note that \\n\\ LO o/Q pL 1 ^ 2 )) 
||no|| L i( R 2) < oo for any t < T* . In Proposition^ in section 2, we prove that if ||no|| L i( R 2) is 
sufficiently small, blow-up does not occur in a finite time. We, however, leave an open question 
whether or not singularity may develop for large L l norm of no. 

Remark 2 Liu and Lor z] 15^ showed global-in-time existence of weak solution to (jl.ip in two 
dimensional case under the assumption (jl.4j) . Since their weak solution satisfies integrability 
n E L 2 (0,T;L 2 (R 2 )) for any T > 0, which is a special case in (jl.7p . their weak solution is, 
in fact, a classical solution if (n ,c ,u ) E H m ~ l (R d ) x H m (R d ) x H m (R d ) for m > 3 as a 
consequence of Theorem{J\ 

The second result is the existence of regular solutions under the assumption that HcoH^oo 
is sufficiently small. 

Theorem 2 Let the assumptions in Theorem[l\hold. We consider the Navier-Stokes equations 
in R 2 and the Stokes system in R 3 in (jl.ip o. There exists a constant 5 > such that if 
ll c o|lz,°°(Rd) < &> then classical solution of (jl.ip exists globally. Furthermore, n and c satisfy 
the following time decay: 

\Wt)\\ Lx(m + \\c(t)\\ Loc(Rd) <C(l + t)-i, d = 2,3. (1.10) 



The proof of Theorem [2] will be given in section 3. Next we study (PP-KSF), namely 

dtn + u • Vn — An = —V • [x(c)nVc] 

d t c + u-Vn = -k(c)n t > 0, x G R d (1.11) 

<9tit + u • Vn — Au + Vp = — rtV0, V • u = 0. 

When the fluid is absent, the Keller-Segel equations with chemical of ODE type, typically 
referred to the angiogenesis system, has been studied in JU0 [6J and [20]: 

d t n- An = -V- [x{c)nVc], d t c = -c m n t>0,xeR d . (1.12) 

In section 4, we show local classical solution of p. lip by the usual iteration method and 
present blow-up criteria of (jl.lip . if a finite time singularity occurs. Now we state the third 
main result. 

Theorem 3 Let m > 3 and d = 2,3. Assume that x(-),K) G C m (R), ||W|| Loo < oo 
for 1 < \l\ < m. Then there exists T* , the maximal existence time, such that if (jiq,cq,uo) 6 
H m (M, d ) xH m+1 (M. d ) xH m (W i ), then there exists a unique classical solution of (jl.lip satisfying 
for any t < T* 

(n,c,u) G C(0,t;H m {R d ) x H m+1 {R d ) x H m (R d )), 

(n,u) £ L 2 (0,t;H m+1 (R d ) x H m+1 ). 

Furthermore, ifT* < oo, then one of the following is true in each case ofR 2 orR 3 , respectively: 

(2D) IMIl2(0,T*;L->(R2)) =°° (!- 13 ) 

3 2 
(3D) IMIz,7(0,T*;L/»(]R3)) + l|n|li2(0,T*;L«'(K3)) = °°» ^ + - < 1, 3 < ^ < oo. (1. L4) 



/3 7 
If fluid equation is the Stokes system for (3D), \\u\\ Lq / t*-lp(r 3 )) * n P-l^P is dropped. 

Last main result is global existence of regular solutions for (jl.lip and their decays in time, 
when ||no|| <j m ^ is sufficiently small. To be more precise, we obtain the following: 

Lj^Z ( 



Theorem 4 Let d = 2, 3 and we consider the Navier-Stokes equations in R 2 and the Stokes 
system in R 3 in (0HU 3 . Suppose that (n ,co,u ) € H m (R d ) x H m+1 (R d ) x H m (R d )) for 
m > 3. T/ien there exists ei = C2(d, ||co||l°°) such that if \\no\\ d < €2, then solutions of 

Li 2 (/x J 

(jl.lip become global and classical. Furthermore, n satisfies the following time decay: 



\n 



L° 



<C(l + t)-\ (1.15) 



This paper is organized as follows. In Section 2, we prove Theorem [T] by obtaining a priori 
estimates. Section 3 is devoted to prove Theorem [2] by adjusting De Giorgi method introduced 
in [20] . We obtain the blow-up criteria in Theorem 3 by using a priori energy estimates in 
Section 4 and The proof of Theorem 4 is presented again by using De Giorgi method in the 
last section. 



2 Blow-up criteria of parabolic system 

We first recall following blow-up criteria for (|1.1|) obtained in [H Theorem 2]: 



(2D) / llVcll^) = oo, (2.1) 



/T* r-T* g 2 

II U II2/3(R3) + / ll Vc lli-(R3) = oo, - + - = 1, 3</?<oo, (2.2) 

where T* is the maximal time of existence. If the fluid equation is the Stokes system, not 
the Navier-Stokes equations, in case of 3D, then the condition on v in (|2.2p is not neces- 
sary and thus it can be dropped. From now on, we denote Lf x = L q (0,T*;L p (M. d )) and 
lT tx = L p (0, T*; L p (M rf )), unless any confusion is to be expected. All generic constants will be 
written by C, which may change from one line to the other and e will be used to indicate some 
sufficiently small positive number. 

Proof of Theorem Q] We argue by contradiction. We suppose that (11 .71) in 2D and (II. 8p in 
3D are finite. We then show that T* cannot be a finite maximal time of existence, which will 
lead to a contradiction. We start with the case of dimension two. 

• (2D case) We first suppose ||n|| L 2 ([o,t*)xR 2 ) ^ s finite. We will show f Q ||Vc||^oodi < oo, 
which is contrary to the blow-up criterion (|2.ip proved in [TJ. In the following, we obtain a 
priori estimates, since our computations are made for any time T with T <T* . We frequently 
use the following type of interpolation inequality 

\\ D /IIlp(r2) < C||/IIl?(]r2)||.d /|| l 7 (R 2), (2.3) 

where < I < k — 1 and I — - = —9- + (1 - 9)(k - 2 ) with < 9 < 1. From maximum 
principle for c and conservation of mass for n, it is immediate that ||c||£°° < HcoHl 00 and 
IM|r°°,i < ||rao||z,i- We note that the convection term, (u • V)c, is estimated as follows: 

t,x 



f T 

(n-V)c|L 2 <C/ IMIl 2 II^ 7 ' u IIl 2 II c IU CX) ||Ac||£2(tt 
*>* Jo 



'0 

<e||Ac|| 2 2 +C e ||u|| 2 r00 ,2||c|||oo ||Vn||2 2 . (2.4) 

t,x L*t x £,& t,x 

Via L 2 -estimate of the heat equation, we have 

\\c t \\ 2 L? +\\Ac\\ 2 L? <C\\cof H i+C\\n\\ 2 L? +\\(u-V)c\\l ? . (2.5) 

t,x t,x t,x t,x 

Combining (12.4=1) and (12.51) together with the hypothesis llnll 2 j2 < oo, we obtain ||Ac|| 1-2 < 00. 
On the other hands, L 2 scalar product for equation of n gives that 

-rIMli 2 + l|Vn||| 2 < C||nVc|| L 2||Vn|| L 2 < -||Vn||| 2 + C||n||| 4 ||Vc|||4 

il/V 4i 

X - 1 11O -^- 1 1 11O 1111O 11O111O 

< - Vn ,2 + Cm L2 Vn L2 c w Ac ,2 < - Vn r2 + Cc L Ac r2 n f2 . 



Using Gronwall's inequality, we have 



Inll 2 

T 



0,2 + ||Vn||? 2 < [Inollia exp ( C||c||f,oo ||Ac||? 2 



Next, testing —Ac with equation of c, we obtain 

Wc\\ 2 r2 + l|Ac|| 2 r2 < / |V-u||Vc| 2 (i2; + C||n||r2 + -||Ac|| 2 



< C||Vu|| L 2||Vc|| 2 4 + C\\n\\ 2 L 2 + ^||Ac|| 2 2 < C||Vu|| 2 2 ||Vc|| 2 2 + C||n|| 2 2 + -||Ac|| 2 2 . 
Again, Gronwall's inequality impiles that 

llVcll 2 ^ + ||Ac||r2 < 00. 

L t,x ^t.x 



The energy estimate of vorticity equation leads to 



1/7 1 

2 j t \H\h + W^Wh < C\\n\\ 2 L2 + - 



Therefore, we obtain Hwlr^ + ||Vo;|r 2 . 2 < 00 via Gronwall's inequality. Finally, we have 

L t,X L t]x 



1 d 



- — ||Ac|| 2 2 + ||VAc||| a < C||VuVc|| 2 2 + C||wV 2 c|| 2 2 + C\\Vcn\\\ 2 + C||V 



2dt 



n\ 



L± 



< CWXluW 2 WXJrW 2 -l-Z^IU/ll 2 IIA/-II 2 4- C \\\1A I 2 llnll 2 4- ^11^7 nil 2 

\ *-^ 1 1 V ti 1 1 ^3 II V C 1 1 ^g — I < — -- 1 1 It 1 1 £00 ||ZAC|| r^ 2 T" <^ || V C|| jj, || 't||^3 T l> || V ll\\ t% 

he bound of ||VAc 
proof for the case ||rt|| L i < 00. 



This gives the bound of ||VAc|| 2 22 and in turn, the bound of II Veil, 2, 00. We complete the 



J t . x 



t,x 

In case that ||n|| t p,i < 00 with - + - < 2 and p > 2, we observe that ||n|| r 2 < 00. Indeed, 



due to interpolation and Holder's inequality, 

T r T p^2 



\n\\ 2 T2 dt<C / ||n||?r 1 ||ra||?; 1 d£<C||n|| , '- : 



Jo x >* 

Hence, with the aid of previous result of L 2 t case, the case | + | < 2 with p > 2 is direct. It 

remains to consider the case ||n|| f p.9 < 00 with — I — = 2 and 1 < p < 2. We note first that 

11 niij. , p 1 



1 d 

2d£"~ ~"l^ 



Mil 2 + l|Vw||| 2 < C||n||ip||u|| r _E l 



2(p-l) 2-p 2p 4(p-l) -, 

<" CWriW t llllll p IIV7l/ll p < CWrtW 3p ~ 2 lllill 3p ~ 2 _L ll\7i/ll 2 
_ ^ 1 1 ' ^ 1 1 LP 1 1 ^ 1 1 r 2 1 1 V U 1 1 , 2 _ ^ 1 1 ' ^ 1 1 r p 1 1 ^ 1 1 r 2 1 — 1 1 ^ ^ 1 1 /~ 2 * 

Therefore, due to Gronwall's inequality, sup ||u|| 2 2 + J Q \\Vu\\ 2 L2 dt < 00. 
Using the mixed norm estimate of the heat equation for c, we have 

||ct|| LP , 9 + IIAcIUp,, < C||c ||^2 + C\\u ■ Vc\\ q LP , q + C||n||| p , g . (2.6) 

Noting that (p ~ — = 1, we compute 







\u ■ Vc\\ q LP dt < 



o 



MIz^pllVcH^pcii < C 



o 



li/ll 9 llrll 2 IIArll 2 dt 
I t*!! r 2p II Mil, 00 II Ml £P 



H§pttt + e/ \\Ac\\ q LP dt<C\\u\\\ x f ||Vn|| L2 Pp " ' dt + e / ||Ac||^,dt. (2.7) 
o Vo L ^ 4 jo jo 

Combining (|2.6p and (|2.7p with sufficiently small e > 0, we have 

\\ru\\ q -4- IIA^II 9 < C 

W^WjjPii ' II llLP'9 — 

Multiplying the equation of n with Inn and integrating it by parts, we obtain 
— fnlnndx + \\VVn\\ 2 L 2 < C\\Vc\\ ^.||Vyn|MI>Ai|l L j i * I 

< CHAclUHIv^ll^ 1 ||Vv^||| 2 < C\\Ac\\l P \\^\\h + \\\^M\h- ( 2 -8) 

Using Gronwall inequality, the estimate (12.8P leads to V-^/re G -L 2 X , which implies ||^||^2 < oo. 
This completes the proof for 2D case. 

• (3D case) Suppose that (|1.8|) is not true. As in 2D case, we then show J Q ||Vc|| it x>(R3) < 
oo, which is contrary to the blow-up criterion (J2.2J) proved in pQ. The proof of the case for 
Stokes system is omitted, since its verification is simpler. We first show Vc G L°°(0,T*;L 2 ) 
and V 2 c G L 2 (0,T*;L 2 ). Testing logn to the equation (fLl"^ . 



— / n log re + 4 
dt 



Vn2 



x(c)VreVc < C 



Vn2 



na |Vc| 



<C 



Vna 



L2 



1 
VA 



llVcll 2 P < e 

L 2 p Lp=t 



Vn2 



+ C e llnllrp II Veil 2p 



(2- 



2p-3 _3_ 

Via Gagliardo-Nirenberg's inequality, we note II Veil 2 P <C||Vc|| r 2 P l|V 2 c|| %, and therefore, 

lyp-=T L, II \\L 

combining (|2.9p . we obtain 

r/ 



rf/ 



n log n + 



V?12 



2 9 2p 

<fllV 2 rll + r \\n\\ 2p -^ WX7 r n2 



L 2 



Multiplying (ll.ip i with —Ac and using integration by parts, 

-— |Vc| 2 + / | Ac| 2 = / uVcAc + / k(c)nAc := I + J. 
We first consider the term J. Following the same computations in (|2.9p - (|2.10p . 

J = - f k'(c)n |Vc| 2 - f kVnVc < e||n3 ||| 2 + e ||V 2 c|| 2 . 2 + C t \\n\\%T* ||Vc|| 
where we used that k' > 0. On the other hand, / is estimated as follows: 

<3-3 /3+3 

/< IMI L/3 ||Vc|| 2^ ||Ac|| L 2 < C||n|| L/3 ||Vc|U ||V 2 c|| L2 3 

L P~ 2 



L 2 > 



(2.10) 



(2.11) 



(2.12) 



< e||V 2 c|| 2 r: , +C f \\u 



L* 



2§_ 
I /3-3 



IVcI 



Summing up (f2TT0|) . (I2J2D and (|2^5]l . 



|(/nlogn+I||Vc||i 2 ) + 



V«2 



L^ 



+ ||V 2 c||* 2 < C e (||n 



L^ 



2p 
2y-3 

LP 



+ « 



20 
0-3 > 



Lrf 
2„ ,- r2 



(2.13) 



||Vc||£ a . (2.14) 



Due to Gronwall inequality, we observe that Vc G L°°(0,T*;L 2 ) and V 2 c G L 2 (0,T*;L 2 ). 
Next we consider the vorticity equation of fluid equations in (|l.ip 



dtu — Au + u ■ Vuj — uj ■ Vu = — V x (nV0). 



Energy estimate shows 



1 d_ 
2~dt 



i II-' 

M| L 2 



+ || Vw|| L 2 = UJ- Vuuj - / V x (raV^)u; := Ki + A" 2 - 



(2.15) 



First, we estimate K^. Following computations as in above, 



K 2 <C |Vn| M < C 



Vn2 



< e 



Vna 



1 
n 2 



IMI 2 P < c 

L 2p L ^x 



Vn2 



L-' 



I 2p-3 _3_ 

In II 2 ||, ,|| 2 p II ry 1 1 2p 

"- TV \\UJ\\ t9 \\vU1 \\ to 

I 1 1 .LP H H_L^ H HL Z 



L^ 



+ e||Va;|| L 2 + C e ||n 



2p 
I 2p-3 
I LP 



6U 



L 2 ! 



On the other hand, for Xi, similarly as in (|2.13p . we show 

2^ 

7v ^ 1 1 r? 1 1 2 1 s~i ll 1 1 jS — 3 11 ||2 

Ki <e\\vuj\\ L 2+C € \\u\\^ \\oj\\ L 2 



(2.16) 



(2.17) 



Adding above estimate together and using Gronwall inequality, we observe that w G L°°(0, T*;L 2 ) 
and VwGL 2 (0,T*;L 2 ). 

Next, considering four cases of | < p < 2, 2 < p < 3, 3 < p < 6, and p > 6 separately, we 
will show that 

||wVc||x,«(o,T*;iP) < C + e||V 2 c|| L9 ( ,T*;Lp)- (2-18) 

The proof of (J2.18P will be given later. Then by the maximal regularity of heat equation, 

II c <IIl<2(0,T*;Lp) + l|Ac|| L 9( 0i T*;LP) < C + e l|V c||i«( ,T*;LP) + C||«||l9(0,T*;LP)> ( 2 -19) 

we obtain ||V 2 c||^g(o i T*;LP) < °°- Now we turn to show that n G L°°(0,T*; Z7) for any r > 1. 
Testing n r ~ l to (|l.ip i and noting \' > 0, we observe that 






N r.r + 



4(r - 1) 



Vn2 



r-1 



L 2 



X I Vcl 2 n r + / xAcn r < C / | Ac| n 



< CIIAcI 



L'< 



n 



Lp- 



p < CIIAcI 



2p-3 _3_ 

r? r ll 2p llr? r ll 2p 



< CIIAcI 



LP 



\n 



2p-3 
r 1 1 2p 



Vn2 



LP II'" IlL 1 

2p 



< a iiAci 



2p-3 
LP 



\n\ 



L r 



+ e 



V«2 



L 2 



(2.20) 



Since V 2 c G L q (0,T*;L p ), via Gronwall inequality, we can prove that n G L oo (0,r*;L r ) for 
any r > 1. Let us choose r > 3 and via (|2.19p we then obtain that || VcH^a^oo < oo. Indeed, we 



note that |[uVc||x t- < C e + e||V 2 c||ir for 3 < r < 6 (see (|2.23p below). Again by the maximal 
regularity of heat equation, we have 

\\ c t\\L 2 (0,T*;L r ) + l|Ac|| L 2 (0jT * ;L r) <C e + e||V c|| L 2( 0jT » ;L r). 

Combined with II Veil r =072 + ||V 2 c||r2 < 00, the above yields to ||Vc|| r 2roo < 00 as desired, 
which is contrary to a blow-up criterion (|2.2p . 

It remains to show the estimate (|2.18|) . By (Vc, u) G LfL^, and (V 2 c, Vw) <E L 2 X , it follows 
that uVc belongs to L 4 (0,T*;L 2 ) n L 2 (0,T*;L 3 ) and so uVc G L q (0,T*;L p ) for' 2 < p < 3. 
The other cases are treated as follows, 
(i) (Case | < p < 2) Setting p* = 3p/(3 — p), we have 

||«Vc|| LP < ||n|| L6 ||Vc|| 6 p ^CIMUIIVcH^llVcll^* <a + e||V 2 c|| L p, 

where we use that 2 < dp/ (6 — p) < p*. Taking L 9 -norm in time variable, it follows that 

uVceLi(0,T*;LP). 

(ii) (Case 3 < p < 6) Using that uj and Vc are in L°°(0,T*; L 2 ), we have 



6-p 3(p-2) 6-p 3(p-2) 2(p-3) 3p 

1 < lU/ii 2p lU/H 2p 11V/-11 <r , iii/ii 2p UmII 2p HVr-H 5p ~ 6 llv 2 fll 5p_<: 
\lp ii 11^11^2 ||«l|^6 II vc IIl°° — ° H"Hl 2 H'HIl 2 ll v Hlx,2 || v HIlp 



6— p 3p 5p— 6 

< rMUiii" 2 ^ - ll\7 2 /~ll 5p ~ 6 < n iii,ii 2(p ~ 3) __u ,= 1 1 xv 2 ^l I 

— ^ II u ll^2 || v HI LP — e II Hl 2 II 1 1 LP 

On the other hands, due to p/(p — 1) < q = 2p/(2p — 3), we note that 



L 

From (1.1)3, we have 



P-2 P 

«ll 2 ,,<llr)ll p ~ :L llT7ll p ~ 1 <r , -l-r , llr)ll' ? f9 9"H 

'* 1 |£2 _^ 1 1 ' ^ II r 1 II "'I I hp — ^ ^ ^ II "Wlp ' \ ) 



j t \W\\h + HVn||| 2 < Ch|| L2 + ||n|| 2 2 . (2.22) 



Combining (|2.21|) and (|2.22|) . we obtain ||w||x,oox,2 + ||Vu|| L 2 < 00. Therefore we obtain 

||nVc|| LP < C 6 + e ||V 2 c|| LP for any t < T* . (2.23) 

(hi) (Case p > 6) We estimate 

p-2 2 2(p-3) 3p p-2 3p 

lliiY7/-ll < IU/II \\\1 r\\ < CIWiiW p \\ii\\ p \\K7r\\ 5p_6 ll\7 2 ,-ll 5p ~ 6 < n\\ii\\ p ll\7 2 ^ll 5p ~ 6 
II u v Hi L p — II II L p II HIl°° — ^ II u IIl°° II II L 2 II II L 2 II HI LP — II IlL 00 II HI L p 

P— 2 p— 2 3p p— 2 3p 

<C||Vu|L 2 2 p ||V 2 u||^ r ||V 2 c|k 5 7 6 < C\\Vu\\ r 2 i ||V 2 c|L 5p T 6 • 

— 11 W L* II II II \\-L p — M u L z II WL*' 

Therefore, suing q = 2p/(2p — 3) and Young's inequality, we have 



||nVc|| LP < C e || Vwllia + e || V^H j*" 6 "* "*' < C £ (l + [| Vw||£ a ) + e ||V 2 c|| LP , (2.24) 

where we used that ( 5 _ 6 -.f 3 _ 4 ^ + < q = 7^3 with p > 6. Therefore, the estimate (|2.18j) is 
also true for p > 6. This completes the proof. rj 

Next we present the proof of existence of regular solutions in case ||wo||.li i s small in 
dimension two. 



Proposition 1 Let d = 2 and initial data, x> k, and <f> satisfy the assumptions in Theorem^ 
Assume further that ||nolnno||Li(R 2 ) + IK^WoIIl^r 2 ) is finite. Then, there exists an e > such 
that if 1 1 wo Hz, 1 < e then the maximal time of existence, T* , is infinite, i.e. T* = oo. 

Proof. Estimates in this proof are a priori, since all computations are made before the maximal 
time of existence, T*. We note that, due to the conservation of mass, sup ||w||^,i < e and 

0<t<T* 

therefore, we have 

||n|| L2 < C\\Vn~\\ L 2 ||Vy^|| L2 < Ce||V^||z 2 - (2.25) 

Multiplying equations of n, c in (|1.1|) with In n, Ac, respectively, we obtain 

± J n tan<b + / 1 vvsf & = - / x( C )A^ < c 4 MM\ WS|». (2.26, 

^H Vc lli 2 + ll Ac ll! 2 < Ce||VV^||L 2 ||A c || L2 + C||Vu|M|Vc|| L 4||Vc|| ia 

< Ce||V^||2 2 + illAcllla + C||Vu|| i2 ||Vc||| 2 + ^||Vu;||£ 2 . (2.27) 

Adding (J2.26P and (J2.27P with the following estimate: 

i a 11 ,|9 I,, , i,9 ^,1, i,9 J-11. , 1 1 o 

-- \\u \\j2 + Vwja < C\\n\\U + - Vw \\U, 
2dt L L L 4 L 



we have 



-71 ( ||Vc|||a + \\oj\\ 2 L 2 + / nlnndx J + ||VVn||£a + ||Ac||£ 2 + ||Vw||| 2 

<C||Vc||i 2 (||V< 2 + ||Vc||| 2 ). 
By Gronwall's inequality, we obtain 

(l|Vc||| 2 + ||w||| 2 + / nlnndx J + / (||W«||z,a + ||Ac|||a + ||Vw||| 2 ) ds 

V J /Jo (2 28) 

< f ll Vc o||| 2 + H^olll 2 + / n lnn dxj exp (C||co||^ 2 ) . 

where we used that ||Vc|| i2 < ||co|| i2 . Next, we estimate Jn|lnn| dx. For simplicity, we set 

D 1 = {x: n{x) < e~^}, D 2 = {x : e~ |a;| < n(x) < 1}. 

A typical argument for dealing with kinetic entropy (see e.g. [7]), we estimate 

/n(lnn)_ = — I nlnn — / nlnn < C I \fn + / (x)n < C I e~~ + / (x)n, 
JD X Jd 2 Jd x Jd 2 J J 

9 1 
where (lnx)_ is a negative part of lnx and (x) = (1 + \x\ ) 2 . We compute 

— / (x)ndx = / nuV(x)dx + / nA(x)dx + / xi c ) n ^ c ^ / ( x )dx. (2.29) 

dt Jm. 2 jr 2 jr 2 Jr 2 
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The term J" R2 nuV(x)dx is estimated by 

/ nu\7(x)dx 
Noting that |V(x)| + |A(x)| < C, we get 



< e||V\/n||L 2 IMlL2- 



nA(x)dx 



+ 



x(c)riVcV{x)dx 



< Ce(l + \\Vy/n\\ r2 ||Vc 



L 2 l|VC|| L 2 y 



Thus, (|2.29j) is estimated as follows: 

/ (x)ndx < Ce(\\Vy/n\\l> + \\u\\ 2 L2 + ||Vc||| 2 + 1). 



d 
dt 



From the it-equation we have 



d_ 
di' 



u 



li2 + ||Vu||| a <Ce\\V^i\\ L 2\\u\\ L 2, 



which gives 



\\ u \\h+ \N u \\ 2 L 2ds <\\u \\ L 2 + Ce ||VV^|| 2 + \\uW\2ds. 
Jo Jo 

We add 2 \fn(lnn)_dx\ to (12T28J) and using ([230]) . (|PT1) we then have 



(2.30) 



(2.31) 



|Vc||^ 2 + IM|| 2 + / n\ \nn\dx + \\u\\ 2 L 2 

+ f QllVVnllia + ||Ac||2 2 + ||Vo;||| 2 + \\Vuf L A ds 
< f ll Vc o|li2 + ||u>o||£a + / nolnn dx) exp (C||c || 2 2 ) 

+C + / (x)n dx + Ce||c ||| 2 + Ce||?/o|lx,2t- 

Therefore, we have L ||V-v/ra||? 2 d£ < C(T), which implies n S l? xt via (|2.25p . Therefore, it is 
direct, due to (II. 7)1 in Theorem [TJ that solutions become regular. This completes the proof, rj 

3 Proof of Theorem [2] 

In this section we present the proof of Theorem [2j We start with the control of L p — norm 
of n under the smallness of ||co[|^oo in next proposition. We use the similar weighted energy 
estimate in [22], which treated the case of x> k are constants and fluid equation is absent. We 
remark that due to the incompressible condition of u, the proof of [221 Lemma 3.1] can be 
applicable to our case and the generalization to non-constant x( c )> K ( c ) is a l so available as long 
as a maximum principle of c holds, i.e. < c < ||co||l°°. 
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Proposition 2 Let the assumptions in Theorem{l\hold andp £ (l,oo). There exists 5 \ = 5\{p) 
such that if \\cq\\ loc < 8\, then n(t) £ L p (R d ) for all t G [0, T*) and 

\\n{t)\\ L v <C = C(p, ||co||z,°°, ||71o||lp), (3.1) 

where T* is the maximal time of existence in Theorem [TJ 

Proof. For a positive 4>(c) such that 4>'(c) > 0, which will be determined later, we obtain 

/ n p 4> = \ n p ~ 1 4> (-u • Vn + An - V • (nxVc)) + - / n p 4>' i-u • Vc + Ac - nk(c)) 

pdt J J p J 

= frf^An- /n p ~VxV-(nVc) + - f n p <j)'Ac 

- J n p ~ 1 (f>x'n\ Vc| 2 - - J n p cp'nk(c) - f n p ~ x ^u • Vn - - I n p cj)'u ■ Vc. (3.2) 

We note that, due to V • u = such that <f)'u ■ Vc = V • (4>u), the last two terms in (|3.2p are 
cancelled, i.e. f n p ~ l 4>u ■ Vn + - f n p <p'u • Vc = 0. Via the integration by parts, we have 

1A /' n P0 + (p_i) f n p - 2 (l)\Vn\ 2 + - fn p (/)"\Vc\ 2 
pat J J p J 

= -2 J n p ~ V Vn ■ Vc + (p - 1) f n p ~ l x<t^n ■ Vc + f n p X <t>'\Vc\ 2 - - f n p (f>'nk. (3.3) 
Noting that the last term in (|3.3|) is non-positive and using Cauchy-Schwartz inequality, 

II fnP^ + l^l / n P-2M V , 2 + I /V/ivd 2 (3.4) 

pdtj 2 J pj 

< ~^-r fn p ^-\Vc\ 2 + (p-l) f n p X 2 (t>\Vc\ 2 + I n P x<t>'\Vc\ 2 . 
We set 0(c) = e^ c ' and we look for <f> satisfying 

-^-r { -^ + (p-l)x 2 <P + X0'<^<p". (3.5) 

p — 1 <p 2p 

Let xi — SU P x( c )- We then see that (|3.5p is satisfied, provided that 

0<C<||co||_£,oo 

(P - l)Xi < ^/3 2 , ||cd|U»Xi < ^-, — ^— r/? 2 ||c |li- < 7^- (3.6) 

op op p — 1 op 

If f} is chosen such that 6p(» — 1)%| = /3 2 and if Xi|I c oI|l°° < 2ip 5 it is straightforward that 
(|3.6p is satisfied. Therefore, if ||co||z,°° is sufficiently small, we obtain 

11 l n P 4>+ V — I ' n p - 2 ( f ) \Vn\ 2 + — f n P 4>"\Vc\ 2 < 0. 
pdt J 2 J 2p J 

Since <fi > 1, it follows that J" Rd n p (t)dx < e^ W c °Wl°° J Rd n p dx. This completes the proof. q 
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We remark that Proposition [2] dose not give control of ||ra||i,°°. Toward the boundedness 
as well as the decay of ||n||i/x>, we modify the approach done in [20], where the degenerate 
Keler-Segel system (|1.12|) was considered. With the aid of incompressibility of the velocity 
vector field u, it turns out that the method of proof in [20] can be adjusted properly to our 
case. Now we are ready to present the proof of Theorem [2J 

Proof of Theorem [2] We first note that solutions are classical, because of Theorem Q] and 
Proposition [2j To obtain a truncated energy inequality for (P-KSNS), we first differentiate 
f(n — K) p + (j) in time variable, where (ft is the function introduced in the proof of Proposition 
[2j Similarly as in (|3.2p . (|3.4p . we have 

~J t J(n- Kf + cf> + (p - 1) j(n - Kf- 2 4>\V{n - K) + \ 2 + -J(n- K) p + <p" \V c\ 2 
= -2 f(n - ijQfpV V(n - K)+ ■ Vc - - f (n - K) 



Y+cftnK 



>o 



+ (p - 1) J(n - Kf-^Vin - K)+ ■ Vc + J (n - K) p +X( p' \V c\ 2 
+ K(p - 1) f(n - Kf~ 2 x4>V{n - K) + -Vc + K I \n - K) p ~ l X <P '|V 'c\ 2 . 
We note that the last two integrands in the equality above are bounded as follows: 
(n - Kf~ 2 ct)xV{n - K) + ■ Vc < — (n - K) p ~ 2 : 0|V (n - K) + \ 2 + 4K(n - Kf~ 2 : 0x 2 |V 'c\ 2 , 

(n - Kf^x^Vc? <(^(n- Kf + + 8k) x^'|Vc| 2 . 



In what follows, we fix p = 2. With the aid of (|3.6p . we have 

~ /(n " K) 2 ^ + \ J d>\V{n - K) + \ 2 + I J (n - K) 2 + ^"\Vc\ 2 
< AK 2 J 0x 2 |Vc| 2 + 8K 2 J ^'|Vc| 2 < SK 2 J (fa 2 + x^')|Vc| 2 . (3.7) 

We set L := sup {&x + X^)- Multiplying the equation c with 16K 2 L, we have 

< C < || CO | | Xy °° 

j f 8K 2 L{c - K)\ + f 16K 2 L\ Vc| 2 < 0. (3.8) 

Summing up (|3.7j) and (J3.8)) . we have 

| (/ |(n - *£* + / 8i^(c - if)?-) +\J^(n- K) + \ 2 

+ \l\n- K) 2 J'\Vc\ 2 + 8iT 2 L j |Vc| 2 < 0. 
Similarly proceeding as in [20], we define 

U (0 = J v(t)J(n-Z V (t))l + J u(t)J(c-^(t)) 2 + :=U 1 (0 + U 2 (0, £ > 0. 
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(3.9) 



Here the auxiliary functions v(t),r}(t), and the range of £ £ [0, 2M] are specified later, where 
u{t) and r/(i) are decreasing in t and M is a fixed number. We then note that 

U'(0 = -2 J v{t)r,(t) f(n-ir 1 {t)) + -2f u(t)r,(t) J \c- £#))+. (3.10) 

Let K\ := sup C??(i)- Repeating similar computations as in (|3.9j) . we observe that 
0<f<2M,t>0 

Jt (/ \ {U ~ ill{t))2+(t) + 8K i L /( c " &#))+) + | / 0|V(n - ^(t))+| 2 
+ ^ /(n - ^(t))^"|Vc| 2 + 8K 2 L | |V(c - £r ? (t)) + | 2 (3.11) 

< -£»/(*) /(" - ^{t)) + 4>dx - ir,'{t) J(c - tr,(t))+. 
We set i^) m i n = min<^>(c) > 1. For simplicity, we define 

E(0 := ^up^ (| i(n - fr(t)& + ^ | l(c - fr(i)&) +\ [ J |V(n - ^(t))+| 2 



1 i r [f c ,^1 y^ ,2 , ^ r 



''mm JO 



/ [(n-Zv(t))W\Vc\ 2 + 8 -^ [ /|V(c -*#))- 

JO J (pmin JO J 



2 



Then, after integrating (|3.1ip in time variable for £ > £o := V X (0) max{||no||.L°°, ||co||z,°°} 5 we 
obtain 

E(0 <-J-(f £r/(t) f{n - Zv(t))+4>dx + / £/(t) lie - &$))+) • (3-12) 

<Pmin \Jo J Jo J J 

Assuming that |?/(t)| < Cv>(t)r](t), which will be confirmed later, we have 

E(®<CZ\U'(Z)\. (3.13) 

We use the Sobolev embedding and then by interpolating we have 

||(n-e#))+|li« +\\(c-Cv(t))+\\h <CE(0, q = 2(d + 2)/d, d>2. (3.14) 

t,x t,x 

For simplicity, we denote 

A:=J Ju(t)$(n-€n(t))+, B:= J f i/(t) = (c - £??(*))+■ 

Interpolating 1 < 2 < g in space and using the Holder's inequality in time, we have 

E7i(fl < A Q ^(e), Da«) < #^(0, « = ^, = j|±|. (3.15) 

Under i/(t)7/(t) > CV(f)« = CV(i) 1+ i, we obtain from (fBTTOl) together with (13TT3J) and (|5151l 
|C/'(0I > CA + CB > CU*E- l a > CC 6/a U«\U'\- e/a , (3.16) 
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where we used that (A + B) a > C a (A a + B a ) for < a < 1. Now we choose the auxiliary 
functions z/ and n by 

*/(*) = (1 + t)" 1 , r/(£) = (l +i )-#. 

By the similar reasoning mentioned in [20], we need to show that U(£) is finite for some £ > 0. 
Indeed, for q = 2(d + 2)/d we have 

U(2£)= / v(t)(n-2£r)(t)) 2 + dxdt+ / i/(t) (c - 2£r)(t)) 2 + dxdt 

JO J{n>2£r)(t)} JO J {n>2£r)(t)} 

j-T i- 2(d+2) 4 

</ K*) / («-^(i))+ d {&}{t)Y*dxdt 

JO J{n>2£r](t)} 

f T f 2{d+2) 4 

+ / v(t) I (c-^(t))+ d (Zr,(t))-*dxdt, 

JO J{n>2£ri(t)} 

where we used that n — £r](t) > £?/(£) on {:r|n — 2£r/(i) > 0}. Therefore, we obtain 

17(20 < / T " (t) 4 fn^d X dt + [ T "® A [c^dxdt 
Jo (^(t))3 7 7o (er?(t))3 7 

fT 2(d+2) fT 2(d+2) 



< C^ 3 ' 



2(d+2) "" I / ll^ll 2(d+2) ' 

d JO d 



(3.17) 



4 2(d+2) 2(d+2) 

<cr 3 (iKii 2 d +iicoii 2 d ; 



where the last inequality in (J3.17P is due to the case of K = in (|3.9p together with (|3.8[) and 
(|3.14p . Therefore, f|3. 17|) implies that U{£) is finite for every £ > as long as U(£) exists. Via 
(|3.16p and (|3.17j) . we observe that 

U'{0 < -CC^U^, e > Co := ^WdlnolUoo, ||co||l~}. 

and it is immediate that U(£) vanish at a finite value £ = A7(£o, ||^o||l 2 , ||co || x, 2 ) - Summing up 
the arguments, we conclude that n(x,t) +c(x,t) < CM(l + t)"3 for t > 0. This completes the 
proof. rj 

Remark 3 The 1? energy inequality (13. 9h is responsible for the time decay rate t~ ' 4 . The 
number coincides to that for the solution of the heat equation with the initial data in L 2 (R. d ). 
We do not know whether or not such decay estimate can be improved, and thus we leave it an 
open question. 

4 Blow up criteria of parabolic-hyperbolic system 

In this section, we consider ([l.lip . which is the case that equation of c is of no diffusion. First 
we construct solutions of (II. lip locally in time in the following class of functions: 

X% := (C([0, T) ; 77 s ) n 7 2 (0, T; 77 s+1 ) x C([0, T) ; H s+1 ) x (C([0, T) ; H s ) n L 2 (0, T; 77 s+1 ). 

Our construction of regular solutions is based on the method of contraction mapping via 
linearizing the equations in an iterative way. Next proposition is the first part of Theorem [3l 
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Proposition 3 Let initial data, \, k, and <f> satisfy the assumptions in Theorem O Then 
there exists T > depending on ||no||#*, ||co||#s+i, ||«o||jf s with integer s > 2 such that a 
unique solution (n, c, u) in Xj, exists. 

Proof. We consider following linearized system, which is defined iteratively (set (n°,c°,u°) = 
(n , c , n )) over R d x (0, T) with d = 2, 3. 



' <W m+1 ) + (n( m ) • V)?i( m+1 ) - An( m+1 ) = -V • [x(c( m ')n (m) VcH], 

5iC (m+l) + ( n (m) . V ) C ( m+1 ) = -/c(c( m ))n( m \ 

^ u (*n+i) + ( U M . V)-u(' m+1 ) - An( m+1 ) + Vp( m+1 ) = n( m )V(/>, 
div - u ( m + 1 ) = 0. 



(4.1) 



• (Uniform boundedness) If the initial data (no, Co, i*o) G -^ s x H s+1 x i? s with integer 
s > [g] + 1, then we show (n^ m >, c^ m ', vS m ') is uniformly bounded in Xj, for some To > 0. Let 
a = (a%, • • • , ad) be a multi-index and \a\ := a± + ■ ■ ■ + a^- Taking D a operator on the first 
equation in (|4.1|) . taking scalar product with D a n ym+l > and summing over \a\ < s, we have 



ii| n (»n+l)||^ + || Vn (m+l)||^ <-Y I D a (u^ • Vn< m+1 ») -D a n( m+ Vdx 
2 at z - / 7 M d 

|q|<s 

+ V / Z/*(x(c( m ))n M Vc (m) ) • L> Q Vn( m+1 W := /1 + J 2 . 
|a|<s 

Using cancellation and calculus inequality, we obtain 

|ii| < 5Z / (# a (u (m) • Vn( m+1 ») - u< m ) • VZ? a n( m+1 )) • L> a n( m+1 >dx 



\a\<s 



< C Y (||Vn (m) || L oo||n (m+1) ||^ + h (m) ||^||Vn (m+1) || L oo)||n( m+1 )||^ < C\\ u (m) ||if.||n ( ' n+1) |&,. 
Using Young's inequality and interpolation inequality, we have 

\h\ < C J2 \\D a {x{c {m) )n^Vc^)\\ 2 L2 + I||Vn( m+1 )||^ 



\a\<s 



< n(\ _|_ ||^( m )|| 2s nln( m )|l 2 ll/-( m )|l 2 _|_ i||Vr7( m+1 )|l 2 

— ( - / l 1 + ll c IliJs+iJII' 1 ll_f/ s ll c ll.f/s+1 + — || V77, \\h b - 



We find that 



dl 



\n 



(m+l),|2 , || V „(m+l)||2 <• /711..M || ||„(»n+l) ||2 

-\-<^{± -\- \\C \\h s + 1 )\\"' ll-ff s ll c ll_H" s + 1- 



Similarly, we have 



dt 



I (m+l)i|2 < (7||V7/( m )||„Jlr (m+1) ll 2 ^ 

|C || H"s+1 _^ ^ || V U H-T2 II C ll^-s+1 



+G(||Vn v ; ||H s + ||n v >\\ L oc)(± + ||c v ; ||# s +i)||c 



(m+l) I 



H a + 1 i 



and 



c//- ' 



. (m+l) || 2 1 ||r7„,(m+l)||2 ^ r<\\-,,( m )\\ \\-,,( m + 1 )\\ 2 _i_ r'll^.M II IU,( m + 1 )ll 
L \\H S 1 N VTi v y \\jjs S (s \\u \\h s \\u ll.ff s + Ml n ll-H" s ll n ll-ff s - 
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Adding above, we have 

Ac||r)( m+1 )|| 2 rr +|U m + 1 )|| 2 J _ 1 + || ? /( m + 1 )|| 2 rr ^4 ||Vr7( m+1 )|l 2 r r 4- llW m+1 ) II?, 
jjAH \\h s i^ \\ c \\h s + 1 ' ll u ll_ff s J ' II v n ll_ff s ^ll vu \\H s 

< C(\\u {m) \\ H s + \\Vu {m) \\ H s + ||Vn (m) || H , + ||n( m )|| LO o + l)x 

ni T7 ( m + 1 )||2 , || (m+l)i|2 , || (m+l)||2 \ , pM , || (m)||2s \||„M||2 || r M||2 

vll''' lli? s T ||C ll_f/ s + 1 ' H U ll_f/ s / T^ ° V 1 T^ ll C llfp+i /ll"' ll-f/ s ll C ll/fs+1 

i /~VllY7,>n( m ) II i ll,>n( m )ll N /i i 1 1 „(m.) 1 1 s \2 , /^r||„(m) ||2 

+G(J|V7r ; ||H* + ||^ ll-L^Jl 1 + ll c lljfs+lj +C-|l n ll# s - 

Gronwall's inequality gives uniform boundedness in Xf- for some To > 0, because 

± dt 
o 



simfllr>( m+1 )|l?, 4- llr(' m+1 )|l 2 4 Il7/(' m+1 ) II?, ^ 4 /" ° II Vr?( m+1 )|l?, 4 HV7/(' m+1 ) II?, 
t'UPvll"' \\H S I 1 1 c — ||j^s+i \ || u \\H a ) < I II v "" \\H a ' II v u ll_ff 

JO 



< (\\no\\h + \\cof Hs+ i + \\u \\h° + CT l J\M 1 / 2 4 M s+2 )) exp (CT 1/2 (1 4 M 1 / 2 ) 
under the hypothesis 

sup (||n( m )||^ 4 \\c {m) \\ 2 H s+i + ||u M |&.) + f ° ||Vn^||?j s + \\Vu^\\ 2 Hs dt < M. 
te[o,T ] Jo 

• (Convergence) To show that {{n^ m \ c> m \ u^ m ')} is a Cauchy sequence in Xj, for some 
< T\ < To, we consider the equations of the difference of solutions 

<9 t (n( m+1 ) - n^) - A(n( m+1 ) - n^) 4 (u (m) • V)(n( m+1 ) - n^) 4 (u (m) - M^-^VnW 

= -V • [x(c< m ')n( m )v c ( m )] 4 V • [x(c(" 1 ~ 1 ))n( m - 1 )Vc( m - 1 )], 
$( c (™+i) _ c (m)) + ( u (m) . v)( c ( m+1 ) - c(™)) 4 (n( m ) - -u^" 1 )) • Vc(™') 

= -/c(c( m ))n( m ) 4 fc(c( TO-1 ))n( m-1 ), 
$( u (m+i) - n( m )) - A(u( m+1 ) - n( m )) 4 (u (m) • V)(u( m+1 ) - ?j( m )) 

4(u (m) - u^" 1 )) • Vu( m ) 4 V(p( m+1 ) - p( m )) = (n(™) - n( m - 1 ))V(/>, 
div (u( m+1 ) -«("*)) = 0. 

Following the arguments similarly in [I], we can prove the convergence. Since its verification 
is rather straightforward, the details are omitted. rj 

To obtain the blow-up criteria in Theorem [3l we derive lengthy a priori estimates. Espe- 
cially, the estimates of ||Vu||£,2£oo is crucial. To obtain a bound of ||Vu||^2 it xj, we first use 
vorticity estimates to obtain L 2 estimates of Vw, and then, we obtain the estimates ||Vw|| r 2,» 

by using the mixed norms Lf% type estimates for Stokes system (see e.g. [Hi]). Then the 
desired blow-up criterion can be obtained by an induction argument. This is the outline of the 
second part of Theorem [3] and now we give the proof. 

Proof of Theorem [3] Since construction of local solution is done in Proposition [3j it remains 
to show the blow-up criteria for (2D) and (3-D). We will show those criteria by obtaining a 
priori estimates as in the below steps for [0,T] for any T <T*, where T* is the maximal time 
of existence. Since f Q \\n\\ 2 Loo dt < oo and ||n(t)||ii = H^ollz, 1 f° r an * £ [0>^1> we n °t e that 
Jo ll n ll|p^ < °° f° r au P < °°- 
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• (Case M 2 ) At first, we consider the case d = 2. 

Step 1-1 (L 2 x H 1 x H l Estimates of (n,c, «)). Testing it to the equation of u, we have 



1 d 

2dt> 



u 



i? 



+ ll^ 7n lli 2 — C|l n llz, 2 ll n llL 2 — ll n lli 2 +C|l' u lli 2 - 



It follows from integration in time that 



sup WKVWl* 

0<t<T 



+ [ \\Vu\\ 2 L2 dt < C(\\uo\\ 2 L 2 + f \\nf L 2dt). (4.2) 

Jo Jo 



'0 Jo 

Consider the equation of the vorticity uj = d\U2 — d^ui. 

d t u + (u ■ V)u; - Aw = V x (nV(j)). 

Multiplying (|4.3p with ui and integrating, we have 



(4.3) 



Id.. ||2 ,,_ no 

--||w|| L2 + ||Va;|| L2 < 
and, therefore, we obtain 



riVcj) x Vuidt 



< C||n|| L2 ||Va;|| L 2 < C||n|| 2 2 + -||Vw|| 2 2 , 



sup ||w(£)|| 2 2 + / ||Vw||| 2 <ii < ||<^o|||2 + / H n lli 2C ^- 

<t<T Jo Jo 



On the other hand, due to mixed norm estimate of Stokes system (see e.g. [Hi]), we note that 
for any p £ (l,oo) 

\ \\Au\\ 2 LP dt < C\\u \\ 2 H 2 + C \\n{t)\\\ v dt + C ||u||| 2 p || Vu||| 2p dt < oo, (4.4) 

Jo Jo Jo 

where we used 



/ 

Jo 



I 11^ II V7 1 1 * A4- s^ MM-" II V7 1 1 ^ 

\U\\ L 2p\\\U\\ L 2pal S ll n lll,oo(o,T;L 2 P)ll VU Wl 2 (0,T;L 2 p)- 



' \\(u ■ V)u\\ 2 LP dt < 
o Jo 

Hence it follows that HViiH^o^L 00 ) < °°- Next, testing n to the equation of n, we obtain 

1 J P i 

— tM*)IIl2 + l|Vn|| 2 - 2 < / x(c)nVcVndx< -\\Vn\\ 2 L 2+C\\n\\ 2 L oo\\Vc\\ 2 L 2. (4.5) 

2 at 7 R 2 4 

Taking V on the equation of c, multiplying Vc and integrating over IR 2 yield that 



ld_ 
2dt 



|Vc||£ 2 < 



VuVcVcdx 



+ 



V(k(c)n)Vcdx 



< CHVtillicollVcHla + C||n|| L o ||Vc|| 2 :2 +C||n|| L oc||Vc|| 2 :2 + -||Vn||| 2 . 
If we add the above two inequalities (|4.5p and (|4.6p , then we have 
d 



(4.6) 



di 



{\\<t)\\h + l|Vc(t)||i 2 ) + ||Vn||£ 2 < CGlVulUoc + C\\n\\U + l)(||n||£ 2 + ||Vc||£ 2 ). (4.7) 
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Using Gronwall's Lemma, we have 



sup (|K*)||£ 2 + \\Vc(t)\\h) + / \\Vn(t)\\i?dt < C(\\no\\h + H Vc oll 
ie(o,T] 



L>> 



x exp(C(r 1/2 ||Vu|| L 2 (0iT;L cx )) + ||n|||a (0)r . Loo j + T)) / (||Vu|| L °o + C||n[|| r «x, + l)dt < oo. 
Step 1-2 (Induction argument) Assuming that for an integer m with 1 < m 

n £ L°°(0, T; H'^ 1 ) D L 2 (0, T; H m ), c £ L°°(0, T : H m ) 
and 



(4.8) 
(4.9) 



u £ L°°(0, T; # m ) n L 2 (0, T; # m+1 ). 
we will show that 

ne L°°(0,T; H m ) C) L 2 (0,T; H m+1 ), c £ L°°(0, T : H m+l ), 

u £ L°°(0, T; H m+1 ) n L 2 (0, T; # m+2 ). 

First, we take D" operator (a = (ai, 0:2) is a multi index satisfying |q| = a\ + «2i | ck| < m-\- 1) 
with the equations of u, scalar product them with D a u and sum over |a| < m + 1, we obtain 

7;-7:\\u\\ 2 Hm +i + \\Vu\\ 2 Hm +i < - Y^ / £> Q ((w • V)u)D a udx + C\\n\\H™\\Vu\\ H m+i. 
z dt 1 1 7^ , 1 ^k 2 

|a|<m+l 

If we use the commutator estimates such that 



H=o- 
then we have 



m+l „ 771+I „ 

J2 D a {{u-V)u)D a u = Y^ [D a {{u-V)u)-{u-V)D a u}D c 

H=(T H=CT 



< CUVnlliooHull^m+i, 



jr||w|||fm+i + ||Vu|| 2 7m+ i < C||Vn||Loo||u|| 2 7m+1 + C||ra||#m. 
Gronwall's inequality gives us that 

u £ L°°(0, T; H m+1 ) n L 2 (0, T; H m+2 ). 

Next, we take D a operator (a = (01,02) is a multi index satisfying \a\ = a.\ + «2j M < tti) 
with the equations of n , scalar product them with D a n and sum over \a\ < m, we obtain 

— rll n lll™ + l|Vn||jjm < - V / D Q ((u- V)n)L> Q rada; + C||x(c)nVc||Hm||Vn||^m. 

cr|<m 

Using integration by parts (choose ctj / 0) and calculus inequality, we have 



V i D a ((u-V)n)D a ndx = V f D a ~ e J ((u -V)n)D a+e 3ndx 



|a|<m ' 



|a|<m ' 
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< C\\u • Vn\\ H m-i \\Vn\\H m < C(\u\l°° \\n\\H m + C||«||wm-i,oo ||Vn|| i 2)||Vn||ff ?i 

< r(\\ii\\ 2 4- ll-j, II 2 Mlnll 2 _l_ IIY7tiII 2 

— ^ III u IIl°° ~t II u \\w m — 1 >°°) \\" \\H m ' 77 v ' t ll_H' m- 

o 



(4.10) 



Also we have 

||x(c)nVc||#m < C||Vc||i4||n||^-m,4 + C||n||L°o||x(c)Vc||.ffm 

< C , ||Vc||| 2 ||Vc||| 1 ||n||| m ||Vn||| m + Ci||n|| L o C (||Vc|| if - +C 2 ), (4.11) 

where C\ and C 2 are absolute constants depending only on ||c||#m, which is bounded in the 
inductive assumption (m — l)-th step. Using ()4.10p and (j4.11|) . we have 

1 ^ llrdl 2 -lIIVtjII 2 < r(\\i,\\ 2 -4-lli,ll 2 'lllnll 2 

o 37 1 1 1 1 H m ' II » n \\H m — < - y \ II U II_L°° ~r ll^llw" 1-1 ' 00 / II ' l \\H m 

+ C(||Vc||| 2 ||n||^ m + ||n||ico)||Vc||fl*» + C7||n|||«, + -||Vn||^ m . (4.12) 

Similarly, taking H m+l scalar product equation of c with D a c and summing over \a\ < m + 1, 

1 d„ 



2dt' 



Hp»+i<- J] / Da (( u -^)c)D a cdx + C\\k(c)n\\ Hm +i\\c\\ H m+i. 

|a|<m+l ' 



Using commutator estimates 



E 

a|<m+l 



D a ((u-V)c)D a cdx 



V i [D a ((u-V)c)-(u-V)D a c]D a cdx 



|o|<m+l 



< C||Vn||Loo||c||^ m+ l + C||Vc||^4||ti||iym + l,4||c||^m + l 



< C\\ Vu\\ L °o \\c\\ 2 Hm +i + C\\ Vc\\la II Vc||^i ||«||^- m+ i II Vn||^ m+1 ||c|| H m+i 



(4.13) 



and Leipniz formula 



||fc(c)n||#m+i < C||fc(c)||#m+i ||n||ioo + C||fc(c)||i,oo ||Vn||#m+i 

< (Ci||c||^m+i +C 2 )||n||Loo +C||Vn||^m+i, 

where C\ and C 2 are absolute constants depending only on ||c||# m bounded in the inductive 
assumption (m — l)-th step, we have 

Tj-jrllclll-m+i - c (\\Vu\\ L ™ + IMIl 00 + 1 )II c IIh™+ 1 + c|l n lll°° 



1 



+ C||Vc|| L 2||u||^m+i||Vu|| Hm +i + -||Vn||^m+i. 



(4.14) 



Adding P~T2]) and (|4TTi|) . we have 

d 



dt 



n \\H m 1 \\^\\H m + 1 ) ' \\^ ^\\ H m + 1 



< C(||Vc||^2||n||| rro + ||n|||oo + ||Vn|| L oo + ||n|| L oo + l)(||n||#m + ||c||# m+ i" 
+C||n|||oo + C||Vc||| 2 ||it||^m+i||Vu||^ m+ i. 
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Since 

1 1 r7 1 1 2 1 1 1 1 2 i 1 1 1 1 2 

I \ / /■» I 11 _1_ iy~i 

\\ vc \\L°°(0,T;L 2 )\\ n \\L 2 (0,T;H m ) "+" W^W L 2 {Q,T;L°°) 
+ ll^ 7n llL 1 (0,T;L 00 ) + \\^ c \\L x (0,T;L 2 )\\ U \\L° o (0,T-,H m + 1 )\\^ U \\L 2 (0,T;H m + 1 ) < °°' 

it follows via Gronwall's inequality that 

\\ n \\L°°(0,T;H m ) + \\ c \\L°°(0,T;H m + 1 ) + l|Vn||x,2(0,T;if*») < °°- 

This completes the proof of 2D case. 

• (Case M 3 ) Next, we consider the case d = 3. 

Step 2-1 (L 2 x i^ 1 x H 1 Estimates of (re, c, u)). Following similar computations as in 2D case, 

we also have the estimate (I4.2p . We recall the equation of the vorticity uj = V x u> 

d t co + (u ■ V)w -Aw = (w- V)u + V x (raV<£). (4.15) 

Multiplying (I4.3D with to and integrating in spatial variables, we have 

1 d 



2dt 



Mli 2 + ll^^lli 2 ^ 



riV<fi x Vtodt 



+ 



uw||Vu;|cfa; 
1 



< C\\n\\ i- 2 1 1 Voj 1 1 f 2 + ||uw||r2||Vu;||r2 < C||n|| r 2 + C||ua;|| r2 + -IIVwII T 2 

— ii ii J-/ ii ii J-/ ii 1 1 J-j 1 1 ii J-/ — ii ii ±j ii 1 1 Jj /l ' ; 



< C||u|| L/9 ||a;|| %b^ +C||n|| L 2 + -||Vu;|| L 2 < C||u|| L/3 ||W|| L2 -t^\\u\\ l2 t -\\ vwy^ 



and, therefore, we have 



sup ||w(t)||£a+/ ||Vw|||adfc< ( ||wo|lia + 
0<i<T Jo V 



2 fi 1 

2 "- 3 IMI? 2 + Cllnll 2 +i||Vw| 



C / ||n||£ a cft J exp (c I 



2/3 



Using the mixed norm estimate of Stokes system, we note that 

r-T i-T i-T 



[ ||Au|| 2 3 cft < CHuolllf! + C I \\n{t)\\ 2 L3 dt + C I 
Jo Jo Jo 



\u\\ L6 \\Vu\\ L edt < oo, 



(4.16) 



where we used 

l-T 



[ ||(u-V)u||iadt< / 

Jo Jo 



u\\ rfillVnll TR dt < \\co 



\D> 



V- 



\l°°{0,T;L 2 )\\^ U} \\l 2 (P,T;L 2 ) 



Again with aid of the estimate of Stokes system, we have 

r-T r-T i-T 



\ \\Au\\ 2 L4 dt < C\\u \\ 2 H 2 + C \\n(t)\\ 2 Li dt + C I \\u\\\ 6 \\V u\\\udt < 
Jo Jo Jo 



(4.17) 



where we used 

-T 



I \\(u ■ V)u\\ 2 L4 dt < / ||«|[|,e||Vt*|||,iactt < ||w||iao( ,r;i5«)II At *lli a (0,r;X,a)- 
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Hence it is direct that HViiH^^TjL 00 ) < °°- Next, testing n to the equation of n as in 2D 
case, we also have (|4.5p . For equation of c, we can obtain (|4.6p without any modification, and 
therefore, it is immediate that (|4.7p . Using Gronwall's Lemma, we obtain 



sup [\\n(t)\\ L2 
te(p,T\ 



|Vc 



||Vn(t)||i 2 ^<C(||n ||i 2 + ||Vco||| 2 ) 
o 



r 



x exp(C(T 1 / 2 ||Vu|| L 2( 0)T;L cx)) + ||n|||2( 0) T ; Lo°) + Z 1 )) / (||Vm||£,°o + C||n|||cx) + l)dt < oo. 

J u 
Step 2-2 (Induction argument) As in 2D case, most of all estimates are the same as those 
given above. Therefore, we just mention different estimates compared to 2D case. Up to 
estimate (|4.12p . all estimates are exactly the same as before and however, the following is 
slightly different form of estimate (compare to (I4.13P ). Indeed, using commutator estimates, 

V / D a ((u-V)c)D a cdx = V I [D a ((u-V)c)-(u-V)D a c]D a cdx 

|a|<m+l \a\<m+l 

< C||Vu||ioo ||c||^- m + l + C||Vc||£3||li||wm+l,6||c||#m + l, 

1 1 

< C||Vm||loo ||c||^ m+ i + C||Vc||£ 2 ||Vc||^i ||Vn||^m+i ||c||#m+i. 
With the above modification, to sum up, we have 

d (\\r>\\ 2 + iwi 2 ^_liivt7|i 2 

-r:\\\n\\H m > \\ c \\H ,n + 1 ) "+" II vn\\ Hm+1 

< C(||Vc|| 2 2 IMItf™ + ll n lll°° + 1| Vti[|Loo + ||n||Loo + 1| VM||^m+i + l)(||n|| 2 ^ m -|-||c|| 2 ^ m+ i)-|-C'||n|| 2 :o c. 
Under the same assumption as (|4.8p and (|4.9p . Gronwall's inequality implies that 



■;?. 



\L°°(Q,T;H m ) + ll c II.L o °(0,T;.ff m + 1 ) + \W n \\L 2 (0,T;H m ) < 



OO. 



This finishes the case of 3D and therefore, proof is completed. 



□ 



5 Proof of Theorem |4] 

In this section, we present the proof of Theorem 01 The following lemma shows weighted 
energy estimate and truncated energy estimate shown in [5] in case that fluid is not coupled. 
It is remarkable that even in the presence of fluid equations, influence of fluid does not appear. 
Indeed, incompressibility causes cancelation of terms involving velocity of fluid, which is a 
crucial observation for the proof of Theorem HI 

Lemma 5 Let </>(•) be an auxiliary function such that (j)'(c) — </>(c)x(c) = and K a positive 
number. Then, the classical solutions to p. lip satisfy the following weighted energy equality 
(j5.ll) and truncated energy equality (|5.2p .- 



d 

~dl 



n 



(c) + 4 



p — 1 
p 



(c) 



V 



n 



(c) 



p/2 
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(p-1) / 4> 2 (c)x(c)k(c) 



n 



p+i 



(5.1) 



dt 



n 



0(c) + 2^ /0(c) 



"TT*. 



p/2 



*-"/(^ 



p+1 



■«" r(c)x(c)fe(c) + (2p - 1)A / ^(c) X (c)fc(c) 



A' 



+ pA 2 / 0(c) 2 X (c)fc(c) 



/? 



A' 



p-i 



(5.2) 



Proof. We note first that 

d n n t (j) — ntp'ct 1 ,_, , _ _ n^'cj 

-- = ^ = -(An- V- (n X Vc)-u-Vn) - —. 

Due to V • u = and </>' = 4>x-, we also observe that 

An - V • (nxVc) - u • Vn = V • (^v(^) - ^-(pu 



and therefore, it follows that 

d n 1 

dt 

Testing p(;|) p_1 to (|5.3p and using the integration by parts, we obtain 



V • I 4>V I - ) - -<t>u 



(5.3) 



d /" /n 



dt 



4>= p 



n 



p-i 



n\ n 



J ■ |0V f - I --</>< 



/?, 



+ - (factdz 



P ( T ) V-(0V(^]]dx-/p 



-j V- f -0uj dx- (p- 1) / f - j 4>xc t dx. 

(5.4) 



We estimate separately each term in ([57 
p-i 



P 



/'(; 



n \ / n 



V-(0V("]]=-p(p-l) /0(5 



p-2 



V|™ 



V • — 0^ dx = p 



> p-i , 
n \ / n 



dx, (5-5) 

7 ) V • (0u) 4 V ( '- ) ■ 4>udx 



, P-i 

n\ _ / n 



p/(2)V«.)* + /v(=)' 



0ttdx = (p — 1) I — ) V • (4>u)dx, 



■n 



(5.6) 



H 



(P _ !) / T 0XQ^ = -(p - 1) 



-(p-1) 



/ I — J 0(x^ • Vc + kn)dx 



^Vvcdx-tp-l)/^) 1 



(fixkndx (5.7) 



p+i 



- j x^a;- 



23 



Adding up (|5.5p - (|5.7p . we obtain (|5.ip . By the Sobolev inequality, for any p with max{l, d/2 ■ 
1} < P < oo, it follows that 



4/ (jv^(c)<(p-l) 



V W 5 



(7(^)^1 



c)C 



7) 



(c) 



L d/2 



(5- 



as long as < c < ||co||z,°° and sup 4> (c)x(c)k(c) := K\ < oo. For the truncated 

0<e<||co||i,oo 

inequality, we proceed similar computations as in those of the weighted equality. We obtain 



dt 



n 



K 



n 
p\ - - K 



p-i 



d ( n 
dt 



ii 



+ 1--K) 41 ct 



p[--k) V- (<t>v( -A -- ± <t>u\ -p f(^-K) n X c t 



n 



1 - K ) 4>XCt- 



We note that integration by parts yields 



\™-k) V ■ ( 4>V [ " ) ) dx = -p{p - 1) 



7? 



K 



p-2 



V 



II 



dx. (5.9) 



With the aid of replacement of j by ( 5 — K ) + K and integration by parts, it is direct that 



-p / I ^ - K) V • \^<t>u) = -pK ! (~-k) V • (<jw) 



o-*J{\ 



K\ V • (dm). 



(5.10) 



It is also straightforward that 



P 



I 



ii 



p-i 



K tixV • (uc) 



[(J--KJ cp x X7.(uc)dx 



> p-i , 
n „ v n 



n 



pjr--K\ r--K + K\ <p X V ■ (uc) - J I - - K\ X V • (uc)dx 
(p-1) f^-K] 4>xV ■ (uc)dx + pK I (^ - K\ 4>xV-(uc)dx (5.11) 



As noticed earlier, due to </>xV • (uc) = u<f)'Vc = V • (4>u), (|5.10p and (|5.1ip are cancelled out 
each other. We also observe that 



77 



pj [■--K\~ 1 n\k = p j '{^-kJ' 1 '(l-K + K^ '-<*•<* 



n 



P [a-* 



P+l 



ii 



X k + 2pK l-r-K) ^ X k + pK z \--K\ <f>* x k, (5.12) 



n 



<P 



p-i 



n 



K I <f>xkn 



\ ~ K ) ^ Xk 
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K j {^-kJ i? X k. (5.13) 



Summing up (|5.9[) - f)5.13j) . we obtain (|5.2p . This completes the proof. rj 

With the help of Lemma 03 the remaining procedures of the proof of Theorem 2] are almost 
identical with those in [20] and, however, we give the sketch of the proof for clarity. 



Sketch of proof of Theorem |4). Let u(t) , r](t) be auxiliary functions, which will be specified 
later and 0(c) given in Lemma 03 We define a truncated energy E(£) by 



E(£):= sup 
0<t<T 



Kc) J+ P Jo 



V 



n 



(c) 



\P/2 



Using (|5.2jl and following similar procedures as in the proof of Theorem [2l it follows that under 
the condition r/(0) > Hnoll-L 00 



,£(£)< sup / 



II 



+ 2 



P -l 
P 7o 



£ri(t)) 0(c) 



V 



■» 



\P/2 

Sv(t) 



'0 



n 



p-l 



e»7(t) <a(c) 



+ (2p-l)ey r?(t) / /(c)x(c)s(c) 



■)? 



^(t) 



+ p£ 2 / r^) 2 / 2 (c) K (c) 

Jo JR d \4>{C) 



II 



tv(t) 



p-l 



By the sobolev embedding it holds that 

p 



n 



£r)(t) 



<CE(H), q = p(d + 2)/d. 



Li([0,T]xR d ) 

On the other hands, we define the level set energy 

r-T 



m - /. J " W \W 



II 



£r)(t) dxdt. 



Differentiating in £, 



u'(0 



T 



pv{t)T}(t) 



II 



p-l 



£rj(t) dxdt. 



Interpolating p — 1 < p < q, we have 



U(0< (J J v{t)^ 

2p 



II 



p-i \ a 

£r](t) I dxdt E{0 6 



for 6 = rf + , a = ^7§r- Let the auxiliary functions v(£),i](£) satisfy that 



p-i 



v{t)^-o) + |t)(£)| < Civ(t)r/(t), 77(c) < C 2 v{t). 



(5.14) 



(5.15) 
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Then E(£), U(£) satisfy the following differential inequalities: 

J U(0 < Cf \U'($)\ a E(S)°, 



c 2 o + c 2 £U(0, 



d 



Firstly we choose v(t) = (1 + t) 1 2 p and rj(t) = (1 + t)" 1 with C\ = (1 + T) 1 2 p. Working 
(|5. 16[) with G(£) = U(£,) a for some < a < 1, we arrive at G(£) vanishing for a finite £i(See 
(step 4) for Theorem 4.1 in [2D] ) under the condition p > ^±2. The same holds for £/(£) and 
we obtain the decay 

||n(t)|U« < C(T)v(t). 

Next, we choose i/(t) = n(£) = (1 + £) with C2 = (1 + T) 2 ? to relax the initial integrability 
of n to p > 2 L +2 I • Up to this point, L°° decay of n(t) depends on T for < t < T. When 
1 1 no 1 1 d is small enough, (15, ip gives 

J_j 2 



n\ 2 (/ , / , , s l n 



< C||noll d . 



from which we have 11^7^11 d+2 < C, and ||n(to)ll ^+2 < C for tn < 1. Now using the 

result for p > 2 (d+2) an< ^ sca ^ e invariance of the norm ||no|| d and ||co||l°° under scaling (JTT9J), 
we conclude ||i-(£)||l°o < j for a uniform constant C. This completes the proof. rj 
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